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Abstract 



Starting from an unintegrated gluon distribution which satisfies a 'unified' equation 
which embodies both BFKL and DGLAP behaviour, we compute the shadowing correc- 
tions to the integrated gluon in the small x domain that will be accessible at the LHC. 
The corrections are calculated via the Kovchegov equation, which incorporates the lead- 
ing ln(l/x) triple- Pomeron vertex, and are approximately resummed using a simple Pade 
technique. We find that the shadowing corrections to xg(x,Q 2 ) are rather small in the 
HERA domain, but lead to a factor of 2 suppression in the region x ~ 10 -6 , Q 2 ~ 4 GeV 2 
accessible to experiments at the LHC. 



Experiments at HERA and the Tevatron have confirmed the rapid increase of the gluon 
distribution as x decreases, which is expected both in the pure DGLAP framework and in the 
BFKL-motivated approach. It is anticipated, however, that at sufficiently small x, this increase 
will be tamed by shadowing corrections. 

The first quantitative studies of gluon shadowing were made by Gribov, Levin and Ryskin [l] 
(GLR) and by Mueller and Qiu 0. It was found that the shadowing contribution, / s had(^ / , k 2 ), 
to the gluon distribution f(Y,k 2 ), unintegrated over the gluon transverse momentum k, is of 
the form 

df shad (Y,k 2 ) n a 2 1 



BY R 2 k 2 



xg(x,k ) , (1) 
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where Y = ln(l/x), nR 2 denotes the transverse area populated by the gluons and g is the 
integrated gluon distribution. The constant C will be specified later. When the shadowing 
term is combined with DGLAP evolution in the double leading log approximation (DLLA) 
then we obtain the GLR equation for the integrated gluon at scale k 2 

dg(x,k 2 ) N c a s , 2 . a| 



dYd\n(k 2 /A 2 ) vr X9{x ' k ] ~ C lVk 2 )J ' ( 2 ) 

where the last quadratic term, which originates from shadowing, is simply the right-hand-side 
of eq. ©. 

The GLR equation effectively resums the 'fan' diagrams generated by the branching of QCD 
Pomerons, which correspond in the GLR approach to gluonic ladders in the DLLA to DGLAP 
evolution. In this approach the triple-Pomeron vertex, which couple the ladders, is computed 
in the leading In k 2 approximation. The GLR equation has stimulated an enormous literature 
[3-18] connected with shadowing effects in deep inelastic and related hard scattering processes. 
One of the important results to emerge from these studies is the computation of the exact 
triple-Pomeron vertex beyond leading lnfc 2 , but staying within the more appropriate leading 
In 1/x approximation. 

The aim of the present study is to take advantage of this precise knowledge of the triple- 
Pomeron vertex in order to perform a quantitative estimate of the gluon shadowing effects 
which can be probed in the low x domain which is accessible at the LHC. To be precise, we 
start from the solution fi of the unshadowed linear equation which embodies both BFKL and 



DGLAP evolution, as well as subleading In 1/x effects p2| . Then we compute the quadratic 
shadowing contribution, —/shad; ho m the solution fi using the more complete triple Pomeron 
vertex. We resum the shadowing contributions using a simple (1,1) Pade-type representation 

f = 1 7 AO. j JT ' ( 3 ) 

1 + /shad/ JL 

and the gluon distribution is then calculated from 

rQ 2 dk 2 

xg{x,Q 2 ) = xg{x,k 2 ) + — f{x,k 2 ). (4) 

The structure of the triple-Pomeron vertex can be extracted from an equation, formulated 
by Kovchegov [p~T| , for the quantity N(r, b,Y). N is closely related to the dipole cross section 
cr(r, Y) describing the interaction of the qq dipole of transverse size r with the proton target. 
To be precise 

a(r,Y) = 2 J d 2 b N(r,b,Y), (5) 

where Y = ln(l/x) and b is the impact parameter for the interaction of the qq dipole with 
the proton. Recall that the dipole cross section is given in terms of the unintegrated gluon 
distribution by |T 

<r,Y) = *^ / ^[l-J (kr)]f(Y,k 2 ). (6) 



In the large Nc limit, the function N satisfies the integral equation [IS 

N(r 01 ,b,Y) = N (r 01 ,b,Y) + F dy 1-2 In T \ N(r 01 , b, y) + 
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which is the unfolded version of eq. (15) of Ref. [17][|. The linear part of this equation cor- 
responds to the BFKL equation in dipole transverse coordinate space. The term containing 
the log denotes the virtual correction responsible for the Reggeization of the gluon, while the 
linear term under the dr 2 integral corresponds to real gluon emission, p is the ultraviolet cut- 
off parameter. The subscripts 01, 02 and 12 enumerate scattering off qq,qg and qg systems 
respectively. The equation resums fan diagrams through the quadratic shadowing term. 

If we rewrite (]?]) in terms of the transformed function 

roo dr 



N(£,b,Y) 



to r 

then the shadowing term has a much simpler form 



J (£r)N(r,b,Y), 



N(£,b,Y) = N {£,b,Y) + 



a s N ( 



o 



71 



Y 



dy 



K 



N{£,b,y) - N 2 (£,b,y) 
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where K is the BFKL kernel in momentum space ||18|| . Here we have made the short-distance 
approximation in which we neglect the terms in comparison to b, so that N is only a function 
of the magnitudes r and b, and N of £ and b. 

We may resum the linear BFKL effects and rearrange (||) in the form 



N{£,b,Y) = N L {£,b,Y) - 



a s N ( 
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(10) 



where Ni is the solution of the linear part of @ with the shadowing term neglected, and G is 
the Green's function of the BFKL kernel 



G(Y — y) = exp 



fasN ( 
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(Y-y)K). 



(11) 



Eq. ( |10"D may be solved by iteration. At large Y(= \nl/x) the dominant region of integration 
is y ~ Y, where G ~ 1, and so the first iteration gives 



N(£,b,Y) = N L (£,b,Y) 
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We now assume that the b dependence can be factored out of N L as a profile function S{b) 

N L (£,b,Y) = S(b)n L (£,Y), (13) 



1 An equation similar to (0) can be found in Ref. juj 



where we use the normalisation 



J d 2 bS(b) = 1. 



Integrating (|TJ) over d 2 b then gives 

h{l,Y) = n L {i,Y) - 

where 
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We use (HP and @ to write ( |l~5"D in terms of the unintegrated gluon distribution. We obtain 



f(Y,k 2 ) = f L (Y,k 2 ) - 



1 



R 2 \ d\nk 2 
where we have used the identities 




fc2 7^ ln [k 2 hiyJ ) 
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/f 7 J o(^) [1 - Mir)} = \ In {^\ Q(i 2 - k 2 



f(Y,k 2 



(18) 
(19) 



Note that the term in square brackets in (|TT| ) is proportional to rii(k,y). Formula ( |i~T| ) is valid 
in the large Nc limit, but for finite Nc we need to multiply the shadowing term by a factor^ 
N 2 ./(Nq — 1) = 9/8. The second term on the right-hand-side of ( |TTD is simply the shadowing 
contribution —f^d °f ©■ Recall that (0) represents the (1,1) Pade approximation of the series 
whose first two terms are given on the right-hand-side of fli~7|) . Moreover, we emphasize again 
that for the linear term fi we use the solution of an equation which embodies both BFKL and 
DGLAP behaviour and which contains major sub-leading effects in logl/x [22]. In Fig. [I] we 
show the results for the integrated gluon xg(x,Q 2 ) obtained from @ and (f|). The shadowing 
term — /jhid i n (0) I s computed from the unintegrated gluon fi of Ref. [p^j , assuming a running 
coupling asik 2 ). 

Several features of the results of Fig. [I] are noteworthy. First we see, as expected, the effect 
of shadowing on xg (x, Q 2 ) decreases with increasing Q 2 . Second, with increasing In (1/x), the 
start of the 'turn-over' towards the saturation limit is evident in the Q 2 = 4 GeV 2 curves. The 
major uncertainty in the predictions arises from the choice of the value of R, as a consequence 

2 Of course the change of the normalisation of the shadowing term by the factor N^/(N^ — 1) does not 
exhaust possible corrections beyond the large Nc limit. One may expect that there will be other contributions 
beyond those leading to equation (Q) . The general structure of those corrections is unfortunately not entirely 
known. One can however expect that after taking them into account it may no longer be possible to obtain 
a closed equation for N since one will have to introduce other independent and more complicated dynamical 
quantities as well. It should also be observed that the non-linear equation (^) does not contain possible effects 
generated by the compound states of more than two reggeised gluons [Ell. 
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of the 1/R 2 dependence of the shadowing term. We have chosen values of R that are consistent 
with the radius of the protonP]. The results of Fig. [I] show that the effects of shadowing are 
rather small and difficult to identify at HERA where, at best, the domain x ~ 10~ 4 — 10~ 3 
at Q 2 ~ 5 GeV 2 can be probed. On the other hand shadowing leads to up to a factor of 2 
suppression of xg (x, Q 2 ) in the Q 2 ~ 5 GeV 2 and x ~ 10~ 6 — 10~ 5 domain accessible to the 
LHC experiments [p3| . 

For completeness we summarize how the Kovchegov equation (^), may be reduced to 
GLR form We start with @ and approximate 1 — Jq(At) by (&r) 2 /4, which is valid provided 
k 2 C 4/r 2 . Then we obtain 

where the integral can be identified with the integrated gluon xg(x,4/r 2 ), where Y = \n(l/x). 
Thus, from ([5]), we have 

/ d 2 bN(r,b,Y) ~ ^r 2 xg{xA/r 2 ). (21) 

Now if (0) is evaluated in the strongly-ordered approximation (tqi 3> Tq 2 ~ r 2 x ) it can be shown, 
using (f|), that it reduces to the GLR form 

d3[X ' Q2) Ncas xg{x,C?) - ^T^i h(x,Q 2 )Y- (22) 



dYd\n(Q 2 /A 2 ) n ttR 2 Q 2 

Comparing with (0) we see that the coefficient^ C — 1. 

In summary, we have quantified the size of the shadowing corrections to xg(x,Q 2 ) using 
a triple-Pomeron vertex which is valid beyond leading InQ 2 , but staying within the leading 
ln(l/x) approximation. The corrections are found to be sizeable for Q 2 ~ 4 GeV 2 and x ~ 
10 -6 — 10 -5 , see Fig. [I]. This domain may be probed at the LHC by observing prompt photon 
production (gq — > ^yq) or Drell-Yan production both at very large rapidities |23| . Of course the 



latter process involves a convolution to allow for the g — > qq transition, which is required for a 
gluon-initiated reaction; consequently somewhat larger values of the gluon x are probed. 
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Figure 1: The effect of shadowing on the integrated gluon distribution xg(x,Q 2 ), at Q 2 = 
4 GeV 2 and 20 GeV 2 . The continuous line is simply the unshadowed xg obtained from fi ||22|| . 
The dashed and dotted lines show the result of shadowing with R = 4 and 5 GeV -1 respectively, 
where the 'radius' R is defined in (|1^) in terms of the proton profile function, S(b). 
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